plotG(vx,vy,char,size,color) ... Hermite Polyline modules

3a73a61f-ea33-42ba-bbbd-b1ca6746a60e 14 4 5 decimal

&[DATE] &[TIME] - &[FILENAME]

vx vy char size color plotG n vx length : plot vx 1 el vy 1 el char size color augment : i 2 n range plot plot vx i el vy i el char size color augment stack : for plot 41 line :

z n RepEnd m z rows : i 1 n range w i el z 1 el : for i n 1 + (m n + (range w i el z i n - el : for i m n + 1 + (m 2 n * + (range w i el z m el : for w 51 line :

data Strip nest stripped in 'R' u data 1 col : v data 2 col : k 0 : i 1 data rows 1 - range v i 1 + el v i el ≤ (v i 1 + el v i el ≥ (¤ R k k 1 + : el u i el v i el 12 mat : 11 line Undefined if 11 line for 21 line Un-nest R v R 1 el : j 2 R rows range v v R j el stack : for v 81 line :

M Unwrap V M 1 col : i 2 M cols range V V M i col stack : for V 31 line :

data sub subspace j 1 data rows 1 - (range Vsub j el i 1 sub range vx i el data j el i sub / data 1 j + el data j el - (* + : for vx 21 line eval : 11 line for v Vsub 1 el : j 2 Vsub rows range v v Vsub j el stack : for v 31 line 21 line :

Round OFF Matrix/Vector M n Round i 1 M rows range m i el M i el n round : for m 21 line :

user original polygon Xo 0.5 - 0.4 - 0.3 - 0.2 - 0.1 - 0.04 - 0 0.04 0.1 0.2 0.3 0.4 0.5 113 mat transpose : Yo 0.007937 0.012 0.022 0.048 0.167 0.556 1 0.556 0.167 0.048 0.022 0.012 0.007937 113 mat transpose : 21 sys

sub 1 :

select subspaces

New Xo sub subspace Yo sub subspace augment : New Xo 1 el Yo 1 el 12 mat New stack : 21 line

isolate New New :

X New 1 col : Y New 2 col : 21 line

<=  populate for more points at will t 01 0.05 range :

Polyline construct

r 2 : X X r RepEnd : Y Y r RepEnd : 31 sys

x c 10 x ≤ (x 1 ≤ (& cases :

Hermite PolyLine t P 16 / t 3^- 3 t 2^* + 3 t * - 1 + (3 t 3^* 6 t 2^* - 4 + (3 t 3^* - 3 t 2^* + 3 t * + 1 + (t 3^41 mat * :

Hermite basis cubics 0...1

x P 1 el x c * x P 2 el x c * x P 3 el x c * x P 4 el x c * 4 1 sys

t j XP t P 1 el X j el * t P 2 el X j 1 + el * + t P 3 el X j 2 + el * + t P 4 el X j 3 + el * + : t j YP t P 1 el Y j el * t P 2 el Y j 1 + el * + t P 3 el Y j 2 + el * + t P 4 el Y j 3 + el * + : 21 sys

U j 1 X rows r 1 + (- (range : ct 1 : j j i 1 t rows range Hx i ct el t i el j XP eval : for ct ct 1 + : 21 line for Hx 41 line :

V j 1 X rows r 1 + (- (range : ct 1 : j j i 1 t rows range Hy i ct el t i el j YP eval : for ct ct 1 + : 21 line for Hy 41 line :

collect the project xy U Unwrap V Unwrap augment : Polyline xy Strip : 21 sys

bSpline modules

u b 0 u 2 - ≤ 16 / 2 u + (3^* 2 - u ≤ (u 1 - ≤ (& 16 / 2 u + (3^* 23 / 1 u + (3^* - 1 - u ≤ (u 0 ≤ (& 23 / 1 u - (3^* - 16 / 2 u - (3^* + 0 u ≤ (u 1 ≤ (& 16 / 2 u - (3^* 1 u ≤ (u 2 ≤ (& 0 cases :

The Gaussian b-spline

Gaussian basis 'b-spline'

2 - x ≤ x 2 ≤ & x b Undefined if

0 b 0.6667

u v Bsp v m el u m - b * (m 1 v rows sum :

X Xo 2 RepEnd : Y Yo 2 RepEnd : 21 line

u ox u X Bsp : u oy u Y Bsp : 21 line

bSpline t 2 X length 1 - (2.1 range : i 1 t rows range vx i el t i el ox eval : vy i el t i el oy eval : 21 line for vx vy augment 31 line :

bx bSpline 1 col : by bSpline 2 col : 21 line

G bx by o 5 black plotG : XoYo Xo Yo + 15 black plotG : 21 line

Xo Yo augment Polyline G XoYo 4 1 sys

Xo 2 - 1.37 - 0.89 - 0.55 - 0.25 - 0 0.22 0.45 0.75 1.3 2111 mat :

Yo 0.519 0.414 - 1.061 - 0.987 - 0.347 - 2.791 1.081 0.765 1.181 0.588 0.466 - 111 mat :

X Xo 2 RepEnd : Y Yo 2 RepEnd : 21 line

Observe:

bSpline is not defined wrt 'X,Y' but it includes 'X'

in the range of 't'. It it sufficient to re-instruct

bSpline to calculate wrt the last XY set.

Both data sets have same Xmin/Xmax, thus

will discretise wrt to the same Xmin/Xmax but wrt

the last XY data set. This genius feature has been

observed in the Gaussian ksmooth.

]]>

u ox u X Bsp : u oy u Y Bsp : 21 line

bx bSpline 1 col : by bSpline 2 col : 21 line

sub 1 :

select subspaces

New Xo sub subspace Yo sub subspace augment : New Xo 1 el Yo 1 el 12 mat New stack : 21 line

isolate New New :

X New 1 col : Y New 2 col : 21 line

<=  populate for more points at will t 01 0.05 range :

Polyline construct

r 2 : X X r RepEnd : Y Y r RepEnd : 31 sys

x c 10 x ≤ (x 1 ≤ (& cases :

Hermite PolyLine t P 16 / t 3^- 3 t 2^* + 3 t * - 1 + (3 t 3^* 6 t 2^* - 4 + (3 t 3^* - 3 t 2^* + 3 t * + 1 + (t 3^41 mat * :

Hermite basis cubics 0...1

x P 1 el x c * x P 2 el x c * x P 3 el x c * x P 4 el x c * 4 1 sys

t j XP t P 1 el X j el * t P 2 el X j 1 + el * + t P 3 el X j 2 + el * + t P 4 el X j 3 + el * + : t j YP t P 1 el Y j el * t P 2 el Y j 1 + el * + t P 3 el Y j 2 + el * + t P 4 el Y j 3 + el * + : 21 sys

U j 1 X rows r 1 + (- (range : ct 1 : j j i 1 t rows range Hx i ct el t i el j XP eval : for ct ct 1 + : 21 line for Hx 41 line :

V j 1 X rows r 1 + (- (range : ct 1 : j j i 1 t rows range Hy i ct el t i el j YP eval : for ct ct 1 + : 21 line for Hy 41 line :

collect the project xy U Unwrap V Unwrap augment : Polyline xy Strip : 21 sys

G bx by o 5 black plotG : XoYo Xo Yo + 15 black plotG : 21 line

search 70 : user bSpline search row : User search user augment : 31 line

Xo Yo augment Polyline G XoYo user 1 1 el user 1 2 el + 50 green 1 5 mat 5 1 sys

User 70 0.1992 1.4206 13 mat